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1. Introduction. The test particle theory of diffusive shock acceleration 
predicts a power law cosmic ray (CR) spectrum in momentum of the form 

N(p) « p“ s , s = -(r+2)/(r-l), ( 1 ) 

where N(p) is the differential number density of CRs at momentum p and s 
is determined solely by the compression ratio, r = u!/u 2 . Here, the shock 
is assumed to be planar, of infinite extent, steady state, discontinuous, 
piece-wise constant, and u 2 (U 2 ) is the flow velocity (measured in the 
shock frame) upstream (downstream) of the shock discontinuity (u, is the 
shock velocity in the lab frame) 2 ’ 3 * 4 * ’ 6 . 

Since strong, classical hydrodynamic shocks have a compression ratio 
of 4, the test particle theory predicts a spectral index of -2, in good 
agreement with the source spectrum of galactic CRs that is inferred from 
observations. This account of the spectral index is widely viewed as 
being superior to previous accounts, which required fine-tuning to get 
the correct spectrum 1 . 

The non-linear theory, however, which allows for the fact that the CR 
pressure may be dynamically significant, accordingly allows the 
compression ratio to be larger than 4. This follows from the facts that 
(a) much of the post shock pressure can be in relativistic CRs, and, for 
a relativistic fluid, the classical compression ratio can be as large as 
7, and (b) a significant amount of energy can be lost from the shock 
front in the form of particles that stream away from the shock before 
being convected downstream of it, and, r may then be arbitrarily large^. 

High production efficiencies (>10%) appear to be as widespread in 
astrophysical particle acceleration as does a spectral index between -2 
and -3, and the possibility that CRs affect the compression ratio needs 
t0 . k e taken seriously. While one could argue that efficiencies of order 
10% are marginally consistent with observation without requiring 
significant non-linear effects, invoking an injection rate that always 
gives an efficiency - 10% in the test particle theory requires a fine- 
tuning assumption about the injection, which is probably objectionable. 

Here we show that the non-linear theory gives rise to spectra that are 
very close to a power law, over most of the range of astophysically 
reasonable shock parameters for supernova blast waves in the hot phase of 
the interstellar medium (HISM), even though r may differ significantly 
from 4. The universality of the spectrum, in this view, is related to the 
nature of the non-linear feedback of the CRs on their own acceleration 
rate, not to the universality of the compression ratio among shocks. 

2. The Basic Calculation . The non-linear, steady-state theory includes an 

explicit description of the waves that couple the particles to the fluid. 

The full set of equations is 


pu = C lt 


( 2 ) 
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pu 2 + P + P + P - Co, 
K e w g *■* 


(3) 

dT dP dP 

w w . e 

- — « u - — + v “ — - RP , 

dx dx ph dx w 


(4) 

RP = C ,T 
w A dx 


(5) 

dP 

T ds = dh - — 
P 


(6) 

3 . . 1 du s „ 3F 3 n 3F 

3 x U s F ^ 3 dx P 3p 3x 3x 3 * 


(7) 


Here, p is fluid density, u is fluid velocity, u g is the velocity of the 
scattering centers, P e is the pressure in energetic particles, P w is the 
wave pressure, P g is the pressure in gas that has not yet been shocked, 
T w is the energy flux due to waves, Vp^ is the phase velocity of the 
scattering waves relative to the background medium, i.e. u— u g , R is the 
damping rate, Tds is the heat per unit mass that enters the fluid via the 
wave damping, h is the enthalpy per unit mass and C^, C 2 » and C 3 are 
constant. Equations (2) and (3) express mass and momentum conservation, 
respectively. Equation (4) is an energy inventory equation for the waves 
with the terms on the right hand side representing, in order of 
appearance, amplification of existing waves by compression, growth due to 
the cyclotron unstable CR gradient, and damping (see 8 and references 
therein). Equations (5) and ( 6 ) describe the energy deposition into the 
thermal fluid by the wave damping, and (7) is the diffusion equation for 
F * 4irp 3 f(dlnp/dlnE), where f(p) is the distribution function and D is 
the diffusion coefficient. These equations can all be solved using a non- 
perturbative, non-separable, analytic technique that proves to be 
extremely accurate in comparison with numerical solutions * • The 
numerical solutions referred to are possible only over a limited 
dynamical energy range. However, the analytic technique exploits the fact 
that in a real situation there is a very large dynamic energy range - ten 
decades for galactic CRs - so it should be even more accurate in the 
latter case than in the comparison with the numerical results. 

We have solved the above eqs. for rapid wave damping (R + °°) and for 
v ^ constant upstream of the shock (Vp^ is set equal to zero downstream). 
Results are shown in Fig. 1 where we have plotted partial pressure = 
(j/3)EF, versus energy, E. Here, j - 1 for relativistic and = 2 for non- 
relativistic particles. Each curve is the partial pressure profile for a 
locus of (Mi, v Dh ) values shown in Fig. 2. Each point on a curve in Fig. 
2 produces essentially the same profile as seen in Fig. 1. The range in 
compression ratios from v ^ = 0 to Ml+ °°, is as follows: (a) 7.1<r<7.5, 
(b) 4. 7<r<5. 3, and (c) 3.9<r<4.6. The dashed lines in Fig. 1 show energy 
spectral indices of - 2 (horizontal) and -2.7, representative of galactic 
CR source and observed spectra, respectively. Our solutions are 
approximately power laws, deviating from the horizontal dashed line by at 
most a factor of 3 over 6 decades of energy. (That they are not perfect 
power laws results from the slowing of the flow upstream of the shock by 
CR pressure). It is clear that a small phase velocity, representative of 
the HISM, allows an extremely wide range in shock strengths, all 
producing source spectra essentially indistinguishable from that required 
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to produce the observed CR spectrum. The uncertainties In the CR source 
spectrum caused by energy dependent leakage exceed the small differences 
between curves (b), (c), and the horizontal dashed line. Even curve (a) 
is marginally consistent with the inferred source spectrum. 

The ratio of Alfvin to Q 

shock velocity is given by, * ' * I 1 I 1 I 1 I 

v A / Ul = 2.2 lOH B/( Ul /n), .. ~ 

where B is the magnetic — 1 -ff\ 

field and n is the number ZD lll\^ — — — 

density. For the H1SM, we $ „ | 

can use B ~ 3pG, n ~ 0.01 {jj ^ iff \ ^ — L_T 

per cm 3 , and shock veloci- OZ. '/ \ c 

ties ~ 1000 km/s. There- _3 |L \ _ 

fore, v A /uj ~ 0.07, and, as I ^ 

seen in Fig. 2, a very wide ^ \ 

range in Mi will produce ^ ~4 ~ \ 

spectra similar to those ■< \ 

seen in Fig. 1. 0- \ 

CD “ 5 \“ 

3 » Scaling Laws * The re- O ^ 

suits of the previous sec- ^ —g — 1 — * L 1 ■ I » I . 1 t I 

tion show that the spectrum 2 4 6 8 101214 

of relativistic particles . rMCDr'N/ / \/\ 

is “quasi -universal" in the L0G ENERGY (eV) 

sense that it depends very ^6* Postshock partial pressure vs. 

weakly on the upstream pa- energy. An energy cutoff = 10l 5 eV has been 

rameters. The reason for used. Each label represents a family of 

this, we argue, is that the curves with a particular ui/u(GeV). 


energy content in relativistic parti- 
cles is a very sensitive function of 
the spectrum, so that their back- 
pressure on the flow regulates their 
spectrum very precisely. This principle 0.10 
is illustrated by the fact that the 
overall compression ratio of the scat- 
terers is always established in a way 
such that the compression ratio of 3 ” 
mildly relativistic particles is vir— ' s 'v 
tually independent of the upstream 
shock parameters. In Fig. 3 the overall 
compression ratio of the shock is plot- 
ted as a function of Mach number for 
the limiting case of zero phase veloc- 0.00 
ity. It is seen that beyond modest Mach 
numbers, r « M 3/4 , regardless of the 
fact that the adiabatic index of the Fig. 

post shock gas (i.e. the relative con- 
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2. Labels same as Fig. 1 . 


tribution of relativistic and non-relativistic particle pressure) varies 
substantially over this range in M x . (The classical result for y = 4/3 
is shown by the dotted line.) This M 3/l * law can be derived 
straightforwardly by merely assuming that the Mach number of the flow 
when it encounters Gev particles, M(GeV), and the compression ratio at 1 
GeV , are ~ constant (we find, for v h = 0 and for 10 < M, < 100 
that 3.3 < M(GeV) < 3.6 and 3.5 < r(GeV) < 3*75). The result then follows 
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from the fact that, because the phase velocity is zero, the upstream 
fluid behaves adiabatically , i.e. (using eq. [2]), 

P « p 5/3 => [M 1 /M(GeV)] 2 - [ Ul /u(GeV)] 8/3 , (8) 

Here u(GeV) is the upstream flow velocity felt by GeV particles 6 . Given 
that M(Gev) and r(GeV) ■ u(Gev)/u 2 are about constant for v p ^ ■ 0, r 
scales as M 3//t+ . 

As is evident from Fig. 3, the shock compression is dramatically 
reduced and stabilized by the introduction of even a small v pltl for the 
scattering waves. For an infinite upstream acoustic Mach number, the 
compression ratio, formally infinite when v p ^ e 0, is reduced to only ~ 7 
when v p ^/ui ~ 0.03. This result can be aerived analytically to high 
accuracy; using the same principle invoked to derive the M 3/4 scaling law 
when v D u vanishes, one can 
l/[2.7 Y tY-DM 2 (GeV)]. 

4. Conclusions . The non-linear 
theory of shock acceleration 
has been generalized to in- 
clude wave dynamics. In the 
limit of rapid wave damping, 
it is found that a finite wave 
velocity tempers the accel- 
eration of high Mach number 
shocks and limits the maximum 
compression ratio even when 
energy loss is important. For 
a given spectrum (see Figs. 1- 
2), the efficiency of relati- 
vistic particle production is 
essentially independent of 
v ph* For tlie three "families" 
shown in Figs. 1-2, the per- 
centage of kinetic energy flux 
going into relativistic parti- 
cles is (a) 72%, (b) 44%, and 
(c) 26% (this includes the 

energy lost at the upper energy cuttoff). Even small v typical of the 
HISM, produce quasi-universal spectra that depend only weakly on the 
acoustic Mach number. These spectra should be close enough to E “ 2 to 
satisfy CR source requirements. 
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